Generating phase shifts from pseudo state energy shifts 
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A simple way to generate low energy phase shifts for elastic scattering using bound-state calcu- 
lations is postulated, validated and applied to the problem of e"'"-Mg scattering. The essence of 
the method is to use the energy shift between a small reference calculation and the largest possible 
calculation of the lowest energy pseudo-state to tune a semi-empirical optical potential. The L = 1 
partial wave for e"'"-Mg scattering is predicted to have a shape resonance at an energy of about 0.13 
eV. The value of Zeff at the center of the resonance is about 1500. 
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One of the most technically demanding problems in 
quantum physics is the scattering problem, i.e. the pre- 
diction of the reaction probabilities when two objects col- 
lide The underlying difficulty lies in the unbounded 
nature of the wave function. This leads to a variety of 
computational and analytic complications that are sim- 
ply absent in bound state calculations, e.g. the Schwartz 
singularities that occur in the Kohn variational method 
for scattering [2, 3]. 

One strategy adopted to solve scattering problems is 
to use bound state methods. There are numerous exam- 
ples of such approaches, one of the most popular being 
the i?-matrix methods that rely on the solutions of the 
Schrodinger equation in a finite sized cavity to determine 
the bejiaviour of the wave function in the interaction re- 
gion [ij. The total wave function is then constructed 
by splicing the inner wave function onto the asymptotic 
wave function. 

However, despite the considerable activity in this area, 
there are a number of problems that are beyond resolu- 
tion. The e"'"-atom problem is a notoriously hard numer- 
ical problem since the atomic electrons tend to localize 
around the positron, thus giving a very slowly conver- 
gent partial wave expansion of the wave function inside 
the interaction region (this should not be confused with 
the partial wave expansion of the asymptotic wave func- 
tion) [iSSHi. For example, the dimensionality of the 
equations to be solved to achieve a given accuracy are 
about 5 times larger for e"'"-H scattering than for e^-H 
scattering. At present, there are a number of positron 
collision problems that are simply inaccessible with ex- 
isting approaches 

This article had its origin in a particular scattering 
problem, namely the determination of the near thresh- 
old phase shifts for positron scattering from the di-valent 
group II and IIB atoms. The dimension of the secular 
equations for bound state calculation on such systems 
are very large, for example a CI calculation of the e+Ca 
state resulted in equations of dimension 874,448 Q. 



Application of the CI-Kohn approach [9*] to determine 
the phase shifts for e^-Mg scattering in the channel 
would result in linear equations that are simply too large 
(w 1,000,000) to be solved by direct methods. Itera- 
tive methods do exist, but there are no robust methods 
that absolutely guarantee convergence [l3|- It is likely 
that the development of an efficient linear solver for the 
class of problems that arise from a basis set treatment of 
quantum scattering would involve a good deal of initial 
effort and experimentation. There is, however, a great 
deal of experience in obtaining the lowest eigenvalues of 
large symmetric matrices [ll| . 

The idea behind the current method lies closest to the 
trivial i?-matrix method 12] which is exploited in Quan- 
tum Monte Carlo (QMC) calculations of scattering [l3| . 
In the QMC, one extracts the phase shift by comparing 
the zero point energy of a finite size cavity to the energy 
of the system wave function in the same cavity. In the 
present method, the phase shift is extracted from the en- 
ergy shift when a reference wave function is enlarged in 
size to account for short and long range correlations. The 
method is applied to e+-Mg scattering in the sym- 
metry and used to predict the existence of a prominent 
shape resonance at 0.13 eV incident energy. This is note- 
worthy since shape resonances are currently unknown in 
e+-atom or e^'^-molecule scattering (l^ . 

Our method proceeds as follows. The initial calcula- 
tion uses a reference CI wave function of product form. 



'i'o = <I>gs(X)0o(r) 



(1) 



The wave function of the parent atom is <I'gs(X) where 
X is the collective set of target coordinates. The wave 
function of the projectile is ^o(r). In general, 0o is a lin- 
ear combination of a finite number of square-integrable 
functions designed to give a good representation of the 
wave function in a bounded interaction region. The en- 
ergy expectation, Eq is given by 



(2) 
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The wave function "^q is then augmented by a very 
large number of additional functions to represent the 
correlations between the projectile and the target con- 
stituents. This augmented trial function is 



'i>gs(X)0o(r) 



TABLE I: Results of CI calculations for the symmetry 
of e+H for a series of J. The number of electron (A'^e) and 
positron {Np) orbitals are listed. The total number of two- 
body functions in the CI basis are in the Nci column. Ener- 
gies are given in Hartree while spin-averaged annihilation rate 
$i(X)0j (r) . (3) (r) are given in units of 10^ s"^ {T for *o is for the tuned 

oo limits 



The trial wave function is used to diagonalize -ffoxact 
giving an energy of Ei . The additional functions do not 
include any that have the same sub-symmetries as those 
comprising ^'q. 

Next, a semi-empirical potential of the form 



ad 
2r4 



1 — exp( 



(4) 



is added to -ffcxact (o^d is the dipole polarizability) . This 
potential only acts on the scattering projectile. Then 
^'o is used to diagonalize TJcxact + Vpoi giving i?poi- The 
parameter p in eq. @ is adjusted until Epoi — Ei. Figure 
[T]is a schematic diagram outlining this procedure. 
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FIG. 1: Schematic diagram showing the strategy used to tune 
the semi-empirical optical potential. 

In the final stage, \E'o is modified to permit (/)o(r) to 
describe continuum solutions, 



continuum 



$gs(X)0, 



continuum 



(r) 



(5) 



The phase shifts of i?oxact + Vpoi are then obtained by 
using \l/continuum as the scattering wave function. 

The method is verified by computing the low energy 
phase shifts and annihilation parameters for s-wave e"*"-]! 
scattering. The reference wave function, Vl/g, consisted of 
the hydrogen atom ground state multiplied by a positron 
basis of 30 £ = Laguerre type orbitals. The energy and 
annihilation rate of vE'o are given in Table [H 

A sequence of successively larger calculations with J 
(the maximum £ value of any orbital included in the 
basis) were done up to J = 12. The energies at a 
given J, {E)j, and annihilation rates, (F)j, are given 
in Table HI A major problem affecting CI calculations 
of positron-atom interactions is the slow convergence of 
the energy with J [1, 0, [3] ■ One way to determine the 
J — > oo energy, (i?)oo, is to make use of an asymptotic 
analysis. It has been shown that successive increments. 



Vpoi)- Also given are the extrapolations to the J 
using eq. (|6} . 
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-0.49772560 


0.00089605 
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0.0040914253 
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7087 


-0.49797276 
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298 
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322 




331 


8239 


-0.49797360 


0.0043795165 
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-0.49797439 


0.005341190 


2-term 


eq. 


m 






-0.49797509 


0.005334089 


3-term 


eq. 


m 






-0.49797509 


0.005264739 



AEj = {E)j — (£')j_ij_to the energy can written as an 



inverse power series 

Ae , Be 



lalislliai, viz 



AEr 



Ce 



(6) 



The J oo limits have been determined by fitting sets 
of {E)j values to asymptotic series with either 1, 2 or 
3 terms. The linear factors, Ae, Be and Ce for the 3- 
term expansion are determined at a particular J from 4 
successive energies ((£')j_3, {E)j_2, {E)j-i and {E)j). 
Once the linear factors have been determined it is trivial 
to sum the series to oo [1, [H, (the J ^ oo limits are 
given in Table |T|). 

The trial function was then used to diagonalize the 
Hamiltonian with an additional polarization potential 
{ad = 4.5 Qq). The energy from this calculation matches 
the 3-term extrapolation in TableUwhen p = 2.0495 oq. 
This value of p is very close to a value of p = 2.051 ao 
that was obtained when a polarization potential of this 
form was tuned to a very accurate phase shift in a semi- 
empirical investigation of e+-H scattering [2]J . The phase 
shifts obtained by integrating the Schrodinger equation 
for the model Hamiltonian with p = 2.0495 ao are de- 
picted in Figure [5] and the level of agreement with the 
close to exact phase shifts could hardly be better. 

Besides obtaining phase shifts, this procedure was used 
to determine the annihilation parameter, Z^s- In this 
case the extrapolation to the J — > oo limits were done 
with an asymptotic series similar as eq. but with the 
leading order starting as Ar/{J + 1/2)^. The ratio be- 
tween the annihilation rates calculated with and ^'i 
can be equated with the enhancement factor, G, for s- 
wave e^-H scattering 2l|. The 2-term extrapolation is 
chosen for the evaluation of the ratio since lack of com- 
pleteness in the finite dimension radial basis will have a 
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FIG. 2: The phase shift for e+-H scattering in the s-wave 
as a function of k (in units of a^^). The sohd hue shows the 
results of the present calculation while the triangles show the 
close to exact phase shifts of Bhatia et al [2^ . 



bigger effect on the 3-term extrapolation [a, ll9|. With 
this choice the enhancement factor becomes G = 5.95, 
which is within 1.5% of the enhancement factor chosen 
by normahzation to a very accurate T-matrix close cou- 
pling calculation 23 1. The predicted Zcs, although 
not shown, lie within 5% of those of Bhatia et al [3] over 
the k e [0,0.7] a^^ range. 

This approach to computing the phase shifts was ap- 
plied to the determination of e+-Mg scattering in the 
Lt = 1 partial wave. The treatment of Mg requires the 
use of a frozen core approximation whose details have 
been discussed elsewhere @, [25| , so only the briefest de- 
scription is given here. The model Hamiltonian is ini- 
tially based on a Hartree-Fock (HF) wave function for the 
Mg ground state. The core orbitals are then frozen. The 
impact of the direct and exchange part of the HF core 
interactions on the active particles are computed without 
approximation. One- and two-body semi-empirical core- 
polarization potentials are then added to the potential. 
The adjustable parameters of the core-polarization po- 
tential are defined by reference to the spectrum of Mg+ 

M- 

The e+Mg CI basis was constructed by letting the two 
electrons and the positron form all the possible configu- 
ration with a total angular momentum of Lt = 1, with 
the two electrons in a spin-singlet state, subject to the 
selection rules, 



max(<'o,^i,^2) < J , 
min(^i,£2) < Ant , 



(7) 
(8) 
(9) 



In these rules £o, £i and £2 are respectively the orbital 
angular momenta of the positron and the two electrons. 
The Hamiltonian for the e+Mg state was diago- 



TABLE 11: Results of CI calculations for the state of 
e^Mg for a series of J (Lint = 3). The threshold for binding 
is -0.83285190 Hartree. Most aspects of the Table are similar 
to those of Table U 
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-0.82525710 


0.029828 


11 
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651006 


-0.82806307 


0.12800208 


12 
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724506 


-0.82817969 


0.14306354 


13 


200 
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798006 


-0.82827695 


0.15662562 


14 


214 


216 


871506 


-0.82835799 


0.16873961 
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1-term 


eq. P 






-0.82871101 


0.338475 


2-term 


eq. ® 






-0.82884022 


0.373490 


3-term 


eq. © 






-0.82886332 


0.315877 



nalized in a CI basis constructed from a large number of 
single particle orbitals, including orbitals up to £ = 14. 
The two electrons were in a spin singlet state. There 
was a minimum of 14 radial basis functions for each i. 
There were 20 ^ = 1 positron orbitals. The largest cal- 
culation was performed with J ~ \A and Li„t = 3. The 
parameter Lint was set to Lint = 3 since this is mainly 
concerned with describing the more quickly converging 
electron-electron correlations [2^ . The secular equations 
were solved with the Davidson algorithm [ll|. 
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FIG. 3: The energy of the "'P" state of e+Mg as a function 
of J. The directly calculated energy is shown as the solid line 
while the J — > 00 limits using eq. ((U with 1, 2 or 3 terms 
are shown as the dashed lines. The Mg -I- dissociation 
threshold is shown as the horizontal line. 

First, it is necessary to get the Mg ground state energy 
in this basis. The limitation Lint = 3 means that only a 
single electron in the model atom can have ^ > 3. Trans- 
lating this to an equivalent CI calculation for the Mg 
ground state resulted in an energy of of £^ = —0.83285190 
Hartree (energy given relative to the Mg^+ core). 
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The energy and annihilation rate of the e+Mg state 
as a function of J are given in Table [TTl Figure [3] shows 
the running estimates of (-E)oo with the J ^ oo extrap- 
olations as a function of J. It is clear that none of calcu- 
lations indicate the existence of a bound state, but the 
energy shift algorithm has to be applied to determine 
whether this is due to the finite basis size. 
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FIG. 4: The elastic scattering cross section for e'''-Mg scatter- 
ing in energy region below the Ps-formation threshold. The 
solid line shows the total cross section while the dashed curve 
shows the partial cross section of the 1=1 partial wave. 

A polarization potential given by eq. ^ with aa — 
72 Oq [21|, |25| (the Mg ground state polarizability) was 
added to original Hamiltonian and p was tuned until an 
energy shift of 0.003588 (= -0.82886332 -I- 0.82525710) 
Hartree was achieved. Figure [3] shows the elastic cross 
section for e+-Mg scattering below the Ps formation 
threshold (at k w 0.25 a^^). The cutoff parameters in 



eq. ([4]) were set to p = 3.032 ao for the s-wave [2l| and 
p = 2.573 flo (derived here) for all the other partial waves. 
The elastic cross section in this energy region is almost 
completely dominated by a p-wave shape resonance with 
its center near k « 0.10 a^^ . 

The value of Z^a at the resonance peak was about 1500. 
This was determined by using an enhancement factor of 
G = 12.5 = 0.3735/0.02983 for valence annihilation. It is 
likely that this is an underestimate since a lack of com- 
pleteness in the finite dimension radial basis usually leads 
to annihilation rates being too small @. 

To summarize, a novel technique has been used to 
demonstrate the existence of a shape resonance in e+Mg 
scattering which has the virtue of being readily de- 
tectable. The phase shift calculations were performed 
using a semi-empirical method j21j| with a tuned poten- 
tial. The tuning of a semi-empirical optical potential to 
features such as bound state energies and resonance po- 
sitions is well known. The unique feature of the present 
approach is that the optical potential is tuned to the 
energy shift of a positive energy pseudo-state. This ap- 
proach to the calculation of phase shifts can be applied to 



other scattering systems which are currently inaccessible 
using existing techniques. 

The calculations upon the e+Mg system were per- 
formed on Linux clusters hosted at the SDSU Compu- 
tational Sciences Research Center and the South Aus- 
tralian Partnership for Advanced Computing (SAPAC) 
and the authors would like to thanks to Grant Ward and 
Dr. James Otto for their assistance. The authors also 
thank Prof. Bob McEachran for a critical reading of the 
manuscript. 
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